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Abstract 

I 

The  linear  theory  for  the  eigenfrequencies  of  the  extended  Pierce  diode  for  resistive,  car 
pacitive,  and  inductive  external  circuits  is  tested  via  computer  simulation  with  the  PDW1 
code,  and  is  verified  to  within  the  accuracy  of  the  simulation  (between  1%  and  .01%  de¬ 
pending  on  the  circumstances).  The  simulations  were  remarkably  inexpensive  in  terms  of 
computer  time,  and  could  have  been  performed  cm  a  personal  computer. 

’/ 

Introduction 


The  classical  Pierce  diode  is  a  theoretical  model  with  relevance  to  several  bounded  plasma 
devices.  It  consists  of  two  electrodes  (either  plates  or  grids)  a  distance  l  apart,  connected  by  a  short 
circuit,  with  a  beam  of  electrons  injected  at  one  electrode  with  constant  velocity  vo  and  absorbed 
at  which  ever  electrode  is  first  encountered.  The  electron  beam  is  neutralized  by  infinitely  massive 
co- moving  ions  (see  Fig.  1).  Pierce  showed  [1]  that  the  behavior  at  this  device  is  governed  solely 
by  the  dimensionless  parameter  a  =  upl/v 0  where  wp  is  the  plasma  frequency  of  the  incoming  beam 
of  electrons.  (Pierce  assumed  stationary  ions,  which  made  no  difference  in  his  model,  but  for  the 
extended  model  presented  here,  the  distinction  is  important.)  Pierce  also  worked  out  the  dispersion 
relation  for  linear  perturbations  about  the  uniform  equilibrium.  This  dispersion  relation  is 


2  9a7 

(*3  +  a3)3 


(l  — 1~*  cos  a)  + 


0(0*  -  Q 

(*3  +  a3)3 


sino-f 


0* 

$3  +  a 3 


=  0 


(1) 


where  0  is  the  normalized  complex  temporal  growth  rate,  t.e.,  $  =  ^(7  +  *w).  Despite  appearances, 
there  is  no  singularity  at  0  =  ±ia.  The  standard  Pierce  model  has  been  studied  by  several  people. 
Pierce  found  the  value  of  a  which  first  becomes  unstable,  Cary  and  Lemons  [2]  correctly  described 
the  dominant  branch  of  the  dispersion  curve,  and  Crystal  and  Kuhn  [3]  have  successfully  verified  the 
linear  behavior  of  the  classical  Pierce  diode  using  the  computer  particle  simulation  code  PDW1  [4]. 
Crystal  and  Kuhn  also  investigated  some  of  the  non-linear  saturation  states  via  the  same  simulation 
methods. 
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Id  an  effort  to  understand  the  effects  of  an  external  circuit  on  bounded  plasma  systems,  Kuhn 

[5]  replaced  the  external  short  circuit  with  a  series  RL  circuit,  and  found  the  linear  modes  of  the 
plasma/circuit  system  (which  I  shall  refer  to  as  the  extended  Pierce  diode).  Kuhn  and  Horhager 

[6]  generalised  this  to  a  series  RLC  circuit,  worked  out  the  dispersion  relation,  and  graphed  the 
eigenvalues  for  the  R,  L,  and  C  elements  taken  separately.  The  modified  dispersion  relation  is 
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(Kuhn  and  Horhager  multiply  (1)  by  (tf1+o1)J  to  clear  the  denominator,  but  this  introduces  spurious 
double  roots  at  0  =  ±ta.) 

It  is  my  purpose  to  use  computer  simulation  (the  PDW1  code)  to  verify  the  linear  theory  of 
Kuhn  and  Horhager  for  specific  values  of  a. 


Choice  of  Parameters 


Because  the  regime  of  interest  is  linear,  and  particles’  deviations  from  their  unperturbed  posi¬ 
tions  is  small,  great  accuracy  can  be  achieved  with  minimal  numbers  of  particles  and  grid  points. 
Most  of  the  simulation  runs  were  performed  with  128  particles  (on  the  average)  and  128  grid  cells, 
with  the  time  step  chosen  such  that  each  particle  moves  roughly  one  grid  cell  per  time  step.  With 
these  modest  parameters,  it  is  possible  to  follow  the  decay  of  a  mode  through  twenty  or  more 
decades  in  amplitude.  The  complete  simulation  parameters  are  shown  in  Table  I.  These  simulation 
parameters  make  the  9  representation  (Equation  2)  most  appropriate  (see  Kuhn  and  Horhager  for 
a  description  of  the  alternate  ij  representation). 

The  simulation  eigenvalues  (values  of  9)  were  obtained  from  the  simulation  as  a  function  of  time, 
using  a  four-time- point  scheme  due  to  Buneman  [7]  when  the  expected  dominant  mode  has  a  single 
complex  frequency,  and  by  a  two-time-point  scheme  when  it  has  a  purely  growing  or  purely  damped 
frequency.  (The  scheme  of  Buneman  can  be  extended  to  yield  the  two  most  dominant  modes  if  both 
are  purely  growing  or  damped;  this  extension  was  used  in  some  cases  in  which  the  modes  were  so 
close  together  in  9  that  the  sub-dominant  mode  did  not  decay  enough  for  an  accurate  measurement 
of  the  damping  rate  of  the  dominant  mode.)  The  eigenvalues  obtained  by  this  scheme  will  be  exact 
only  when  there  are  no  competing  modes.  When  competing  modes  are  present,  the  error  in  the 
computed  values  of  the  eigenvalues  will  be  proportional  to  the  amplitude  of  the  competing  mode. 


In  order  to  estimate  the  precision  of  the  simulation  results,  the  eigenvalues  given  by  the  sim¬ 
ulation  (estimates  of  $)  are  plotted  as  a  function  of  time  and  as  a  table  of  numbers.  Typically, 
these  plots  begin  by  fluctuating  strongly,  then  settle  down  to  a  fairly  constant  value,  then  begin  to 
fluctuate  again,  either  due  to  non-linearity  (when  the  mode  is  growing),  or  round-off  noise  (when 
the  mode  is  damped).  Figure  2  shows  an  example  of  this  behavior  for  the  classical  Pierce  diode  at 
a  =  t/2  (the  first  entry  in  Table  II).  The  minimum  fluctuation  about  this  constant  value  is  given 
in  the  tables  of  results  as  an  estimate  of  the  precision  of  the  simulation. 

The  errors  listed  in  the  tables  are  not  meant  to  reflect  the  accsracy  of  the  simulations,  but  rather 
the  precision  with  which  the  simulation  produces  its  result.  The  actual  error  also  includes  systematic 
errors  which  can  be  only  roughly  estimated  knowing  the  order  of  accuracy  of  the  algorithms.  The 
particle  integrator  is  accurate  to  second  order  in  all  cases,  but  the  circuit  advancement  algorithm  is 
of  different  orders  for  different  external  circuit  elements.  These  circuit  advancement  algorithms  will 
be  discussed  as  the  results  for  the  respective -circuits  are  discussed. 

The  a  values  t/2,  3t/2,  and  5 t/2  were  chosen  for  simulation.  These  values  of  a  have  nearly 
maximal  separation  between  the  dominant  mode  and  the  next  most  dominant  one,  thus  ensuring 
that,  in  most  cases,  after  a  reasonable  amount  of  time,  only  the  dominant  mode  will  be  significant. 
These  values  for  a  have  the  additional  advantage  that  they  sample  each  of  the  three  different  types 
of  behavior  of  the  classical  Pierce  diode,  namely  dominant  modes  which  are  monotonically  damped, 
monotonically  growing,  and  oscillatory  growing. 

Initially,  the  beam  fills  the  simulation  region  uniformly  except  for  a  sinusoidal  perturbation  in 
position.  This  perturbation  is  chosen  to  give  the  dominant  mode  an  advantage  over  the  other  modes. 
The  density  profile  for  the  dominant  mode  is  never  truly  sinusoidal,  but  the  modes  are  similar  to 
sinusoids. 

The  external  circuits  also  require  initial  conditions.  Recall  that  the  average  current  of  the 
beam  (which  is  quite  large)  is  neutralised  by  the  co- moving  ions  in  order  to  prevent  a  large  offset 
voltage  in  the  resistive  case,  and  the  build-up  of  charge  in  the  capacitive  case.  Thus  the  external 
circuit  must  respond  only  to  the  perturbation.  The  initial  conditions  for  the  resistive  and  inductive 
circuits  present  no  special  problems,  since  only  the  initial  voltage  and  current  need  be  specified, 
and  regardless  of  how  they  are  chosen,  the  dominant  mode  will  eventually  take  over.  The  initial 
conditions  for  the  capacitive  case,  however,  have  an  additional  complication.  The  charge  on  the 
emitting  electrode  minus  the  charge  Q  on  the  capacitor  is  a  constant;  hence,  if  the  unperturbed 
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|  state  (uniform  beam,  sero  voltage  across  the  diode)  is  to  be  an  allowed  state  of  the  system,  then 

|  the  initial  charge  on  the  capacitor  must  be  chosen  to  be  equal  to  the  initial  charge  on  the  emitting 

electrode.  This  condition  can  be  satisfied  by  choosing  the  initial  condition  for  the  charge  on  the 
external  capacitor  consistently  with  the  initial  electric  field  at  the  emitting  electrode  (which  is 
proportional  to  the  charge  on  the  electrode  according  to  Gauss’  law).  This  choice  is  essential,  but 
fortunately  requires  only  a  little  algebra. 

Simultion  Results 

(a)  Capacitive  Case 

When  the  external  circuit  consists  of  a  capacitor,  it  is  to  be  expected  that  the  smaller  the 
capacitance,  the  more  damped  each  of  the  inodes  is,  since  the  limit  C  — *  0  reduces  to  the  open 
circuit  case  which  shows  no  instability.  The  stabilising  effect  of  the  external  capacitor  can  be 
seen  in  Fig.  3  (taken  from  Kuhn  and  Horhager  —  who  denote  the  normalised  circuit  quantities  with 
overbars,  whereas  here  the  overbars  are  omitted).  Figure  3a  shows  the  real  parts  of  the  growth  rates 
as  a  function  of  a  for  several  values  of  1/C,  and  Figure  3b  shows  the  imaginary  part  of  the  same 
growth  rates. 

The  simulation  results  for  the  capacitive  case  (including  the  C  *  oo  case,  which  is  the  classical 
short  circuit  Pierce  diode),  are  shown  in  Table  II.  The  systematic  error  is  of  second  order  in  the  time 
step  (roughly  0.01  percent),  since  the  circuit  is  an  instantaneous  boundary  condition  on  the  solution 
of  Poisson’s  equation.  Confirmation  of  the  linear  theory  is  plain. 

(b)  Resistive  Case 

One  might  expect  an  external  resistance  to  have  the  seme  sort  of  stabilising  effect  that  a 
capacitor  has,  but  this  is  not  the  case.  Figure  4  (also  from  Kuhn  and  Horhager)  shows  the  dispersion 
plots  for  several  resistance  values.  An  external  resistor  does  not  reduce  the  real  parts  of  «//  growth 
rates,  but  rather,  for  large  resistances,  tends  to  reduce  their  magnitudes,  bringing  them  closer  to 
sero,  which  implies  the  existence  of  a  static  mode  in  the  limit  R  —  oo.  This  seems  to  be  a  paradox, 
since  the  R  —  oo  limit  should  be  an  open  circuit,  just  as  with  the  C  — »  0  limit.  The  reason  for  this 
apparent  paradox  is  the  extra  initial  condition  which  was  imposed  on  the  capacitive  case  in  order  to 
make  the  unperturbed  state  accessible.  For  large  values  of  the  resistor,  the  circuit  becomes  an  RC 
circuit  perturbed  by  the  plasma  where  the  capacitance  comes  from  the  parallel  plate  capacitance  of 


the  confining  electrodes.  The  relaxation  time  of  this  RC  circuit  grows  in  proportion  to  A,  so  the 
growth  rates  (or  at  least  the  dominant  growth  rate)  must  approach  zero  as  1/A. 


The  results  for  the  resistive  case  are  shown  in  Table  III.  Again,  the  confirmation  is  clear, 
although  the  discrepancy  between  theory  and  simulation  is  about  one  percent.  The  method  used 
for  advancing  a  resistive  circuit  in  time  is,  however,  only  first-order  accurate  in  the  time  step,  so  an 
error  of  a  percent  is  to  be  expected. 

(c)  Inductive  case 


An  external  inductor  is  also  a  stabilising  influence  in  most  cases.  As  one  might  expect,  the 
higher  frequency  modes  are  affected  the  most.  Figure  5  (again  from  Kuhn  and  Horhager)  shows  the 
dispersion  plots  for  several  inductance  values. 

There  is  one  exception  to  the  stabiizing  influence  of  the  external  inductor.  Surprisingly,  it  occurs 
in  the  most  stable  regime  of  the  classical  Pierce  diode  —  the  a  <  x  regime.  A  new  mode  is  observed, 
which  is  stable  for  some  values  of  L,  and  unstable  for  others.  Figure  6  shows  the  real  and  imaginary 
parts  of  the  growth  rate  of  this  mode  for  a  =  r/2  versus  the  logarithm  of  L.  The  imaginary  part 
of  the  growth  rate  reveals  that  this  mode  is  in  essence  an  LC  oscillation  perturbed  by  the  plasma, 
where  the  capacitance  comes  from  the  vacuum  capacitance  of  the  bounding  electrodes.  As  the  value 
of  L  decreases,  the  imaginary  part  of  the  growth  rate  (s.e.,  the  frequency)  increases  as  the  inverse 
square  root  of  L.  Expanding  the  growth  rate  Ae  0  for  small  a  and  small  Ae  6  yields 
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Since  only  the  magnitude  of  the  growth  or  decay  depends  on  a,  and  not  whether  the  mode  is  stable 
or  not,  a  transit  time  effect  is  indicated  rather  than  a  true  plasma  interaction. 

The  principal  results  for  the  inductive  case  are  shown  in  Table  IV.  Once  again,  the  theory  is 
confirmed  to  high  accuracy. 

As  mentioned  above,  the  LC  mode  has  the  interesting  property  that  it  can  be  unstable  in  the 
a  <  w  regime,  which  is  stable  for  all  other  external  circuits.  Note  that  Fig.  6  shows  the  beginning  of 
an  infinite  progression  of  stable  and  unstable  regions  as  L  — •  0.  These  modes  are  somewhat  difficult 
to  simulate,  since  the  imaginary  part  of  0  (the  frequency)  varies  as  l/\/T,  and  so  becomes  very  large 
for  small  values  of  L.  The  difficulties  can  be  overcome  by  a  very  small  time  step,  however,  and 


Table  V  shows  the  results  for  eight  values  of  L  (each  at  a  different  extremum  of  the  real  part  of  the 
growth  rate  curve  shown  in  Fig.  6).  The  agreement  is  surprisingly  good. 

Summary 

The  extended  Pierce  diode  of  Kuhn  and  Horhager  has  been  simulated  with  great  accuracy, 
confirming  their  linear  theory  in  all  cases. 
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System  length 

1 

Number  of  grid  cell* 

128 

Time  step 

1/128 

Number  of  time  steps 

1024 

<0 

1 

qjmt 

-1 

mi 

00 

*0 

1 

Injected  electron  current 

—a* 

Background  current 

a* 

Injected  electron  flux 

128 

Table  I.  Simulation  parameter* 
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CAPACITIVE  CASE 


c 

Theory 

Simulation 

a 

oo 

-1.9269 

x/2 

20 

-2.0120 

-2.0122  ±  .0002 

10 

-2.0926 

—2.0928  ±  .0002 

5 

-2.2420 

-2.2423  ±  .0003 

oo 

0.76575 

0.76591  ±  .00001 

3t/2 

20 

0.57219 

0.57238  ±.00001 

10 

0.38749 

0.38770+ .00001 

5 

oo 

0.5892 
+  1.4289* 

0.5898  ±  .0002 
+  (1.4296  +  .0002)* 

5*/2 

0.4602 
+  1.7399* 

0.4605  ±.0001 
+  (1.7404  ±.0002)* 

10 

0.3389 
+  1.9831* 

0.3393  ±  .0002 
+  (1.9833 +  .0001)* 

5 

0.1162 
+  2.3540* 

0.1167 +.0001 
+  (2.3543 +  .0001)* 

Table  II.  Values  of  the  complex  growth  rate  9  for  different  values  of  a  and  C. 


3 


RESISTIVE  CASE 


R 

Theory 

Simulation 

a 

0.05 

-1.7642 

-1.7633 

*/2 

0.2 

-1.3587 

-1.3564 

0.5 

-0.8746 

-0.8728 

2.0 

-0.2908 

-0.2905 

0.05 

0.6407 

0.6417 

3x/2 

0.2 

0.4348 

0.4353 

0.5 

0.2674 

0.2677 

2.0 

0.09229 

0.09235 

0.05 

0.1173 
+ 1 .3827s 

0.1 170  ±.0005 
+  (1.386  ±. 001  )i 

5*/2 

0.2 

-0.6762 
+  0.9299« 

-0.6810  ±  .0005 
+  (0.931  ±.001)i 

0.5 

-0.33648 

-0.33640  ±  .00005 

2.0 

-0.06701 

-0.06778  ±  .00008 

Table  III.  Values  of  the  complex  growth  rate  0  for  different  values  of  a  and  R. 
Errors  not  shown  are  smaller  than  accuracy  given. 
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INDUCTIVE  CASE 


L 

Theory 

Simulation 

o 

0.01 

-0  0775 
+ 10. 0537. 

-0.0777 

10.0563. 

0.1 

-0.56102 
+  3.25597i 

-0  56109 
+  3  25605. 

»/2 

1.0 

-0.08775 
+  0.81437. 

-0.08774 
+  0  81438. 

40.0 

-0.0021749 
+  0.12622. 

-0  0021747 
+  0.12622. 

0.01 

0.74385 

0.744  ±  001 

0.1 

0.61736 

0.61748 

3w/2 

1.0 

0.33751 

0.33758 

40.0 

0.06918 

0.06919 

0.01 

0.51971 
+  1.2920. 

0.520  ±  .001 
+  (1.2925  ±  0005)i 

0.1 

0.25327 
+  0.8424. 

0.25346 
+  0.8427. 

5s/2 

1.0 

0.04338 
+  0.3419. 

0.04342 
+  0.3420. 

40.0 

0.001183 
+  0.05635. 

0.001184  ±  .000002 
+  0.05638. 

Table  IV.  Values  of  the  complex  growth  rate  0  for  different  values  of  a  and  L 
Errors  not  shown  are  smaller  than  accuracy  given. 


LC  MODE 


L 

Theory 

Simulation 

0.00112 

-0.02878 
+  29.92419i 

-0.02996 
+  29.92480i 

0.02521 
+  26.77038» 

0.02536 
+  26.77048i 

H 

-0.03906 

+23.62486i 

-0.03921 
+  23.62513i 

wm 

0.03147 
+  20.47098i 

0.03159 
+  20.470981 

0.0034 

-0.05696 
+  17.21925t 

-0.05709 
+  17.21943i 

■■ 

0.04128 
+  14.22469i 

0.04136 
+  14.22464i 

wm 

-0.10147 
+  10.95935t 

-0.10156 
+  10.95951i 

wm 

0.05809 
+  7.831901 

0.05815 
+  7.83185i 

Table  V.  Complex  growth  rate  for  a  —  r/2  and  various  small  values  of  L. 
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